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Recent first principles calculations of the phonon dispersion curves in MgB2 have suggested the 
presence of anomalies in some of the curves, particularly in the longitudinal acoustical (LA) branch 
in the F to A direction. Similar behavior has been observed in numerous other superconductors 
with F c 's higher than those of standard electron-phonon BCS superconductors. Phenomenological 
calculations of the r — » A LA dispersion based on both an acoustical plasmon and a "resonant 
polarization" mechanism are given here to emphasize the importance of these similarities. 
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Single crystals of the new superconductor MgB 2 large 
enough for inelastic neutron scattering experiments have 
not yet become available, but first principles calculcu 
tions of the phonon dispersion curves have appeared BEl 
In Ref. ^, a well-defined anomaly appears in the dis- 
persion of the LA mode in the r to A direction. This 
anomalous behavior- is remarkably similar to that found 
in transition metalEI and transition metal carbideo su- 
perconductors with high transition temperatures. In- 
deed, it has often been considered a signature of high- 
temperature electron-phonon superconductivity in these 
systems. Here, we call attention to and emphasize the 
potential importance of this similarity by carrying out 
model calculations in which the anomalies are produced 
both by an acoustical plasmon (AP) mechanism and by 
a parameterized, q-dependent, interband contribution to 
the dielectric response function. 

Many superconductors with even moderately high 
transition temperatures show some type of anomaly in 
the phonon dispersion curves when compared to their 
low-T c or nonsuperconducting analogues. In this connect, 
tion, the lattice dynamics in high-T c cuprate materialsu 
is currently being revisited. It now appears that there 
are effects in the optical modest] and in the electronic 
structured that correlate with hole doping and that may 
be related to superconductivity. 

We are mainly interested in anomalies in the LA dis- 
persion curves because they are distinctive and preva- 
lent, occurring in systems that have no optical modes as 
well as those that do. They were first observed and dis-i 
cussed for Nb (T c ~ 10K) more than thirty years agdjJ 
and somewhat later in the transition metal carbidesQ. 
Figure [j] shows the LA dispersion curves for NbC, TaC, 
and HfC, which have transition temperatures of approxi- 
mately 11K, 11K, and <1K, respectively. The anomalies 
are obvious for the first two materials but do not occur in 
HfC. Many other examples of this type of behavior could 
be given, notably the elemental materials Nb, V, and Ta 
where only acoustical modes occur. Also shown on Fig. 
|l| are the results from Ref. § of recent first-principles cal- 
culations for MgB 2 (T c = 38. 5K). While only four points 
were calculated in the r to A direction, they clearly show 
that the behavior is anomalous in the sense being used 
here. 



Nakagawa and WoodsH found they could reproduce 
their phonon data for Nb with a Born- von Karman model 
utilizing numerous force constants. The need for such a 
complex model was thought to be indicative of a strong 
electron-phonon interaction. Ganguly and Woods ob- 
served that the anomalies may be produced by the inter- 
action of acoustical plasmons and phonons. While it is 
difficult to reconcile this suggestion with the electronic 
band structure of Nb, as pointed out by Ruvaldslll, the 
notion that there is a connection between AP-like elec- 
tronic modes and phonons persisted. The resonant elec- 
tronic polarization arising from the double shell model of 
Weber et al£3 is somewhat related to the acoustical plas- 
mon mechanism in that an additional mode which inter- 
acts with the "baae^i phonons is introduced. First princi- 
ples calculationsOO that began to appear subsequently 
could replicate the anomalies but the assignments to spe- 
cific mechanisms were often in conflict. However, the im- 
portant roles of electronic polarization and details of the 
band structure were clear. 

The acoustical plasmon concept was introduced by 
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FIG. 1: Phonon dispersion curves along (q, 0, 0) for NbC, 
TaC, and HfC; the anomalous behavior in the first two is 
evident. Data from the recent first-principles calculations for 
MgB 2 (Ref. |) are also shown. Lines are guides to the eye, 
r.l.u. means reciprocal lattice units. 
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PinesQ in 1958 and subsequently, .acplied to supercon- 
ductivity by a number of authoralilllj'E£(EZl. It is thought 
that APs may occur in systems in which heavy and light 
mass carriers coexist and the light carriers screen the 
heavy ones to produce an acoustical (u — > as q — > 0) 
excitation. Transition metal systems with the d electrons 
playing the role of heavy particles and the s electrons 
the light ones have been long seen as good candidates. 
In other systems, holes may be the heavy particles and 
electrons the light ones. In the simplest acoustical plas- 
mon treatment of superconductivity in transition metals, 
the plasmons either replace entirely or simply comple- 
ment the phonons in a standard BCS-type theory. An ex- 
tensive criticism of AP-mediated superconductivity has 
been given based,on system stability and other restrictive 
requirements! 1 tI 18 !. Nevertheless, while to our knowledge 
APs have not been shown to exist in any bulk mate- 
rial (see also Ref. ^|), the concept is an intriguing one 
and based on sound theoretical arguments regardless of 
any role it may have in superconductivity!. We note that 
acoustical plasmons have been proposedEj as existing in 
MgB2, although not in the context of the present lat- 
tice dynamical discussion and in a much higher frequency 
range. 

Magnesium diboride has the AIB2 structure consisting 
of hexagonal layers of Mg intersper sed bet ween graphite- 
like B layers. Several calculationsEilE3£3 establish that 
the occupied bands consist of a bands formed from boron 
sp 2 hybridized orbitals and ir bands formed from the 2p z 
orbitals directed along the c axis. Interactions with the 
Mg 3s orbitals result in the a bands not being fully occu- 
pied, as in graphite, which creates holes at the top of the 
bands. The various calculations are in good agreement in 
showing a weakly dispersing band in the r to A direction 
lying just above the Fermi surface. Voelker et aZ.,E3 have 
provided analytical fits to the bands they calculated, with 
the holes in the a bands having high effective mass for 
motion along the k z direction. This led those authors to 
investigate the possible existence of acoustical plasmons. 
They found evidence that they do exist but the validity 
of their calculation was subsequently questioned by Ku et 
al.t3 We emphasize that in any case the APs considered 
in these calculations are in a far different frequency range 
than those needed here, as will be seen in the following. 

We write the q- and w-dependent dielectric function 

as 
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and look for its zeros in a standard approach, a 1 and a h 
are the polarizabilities of the two distinct groups of light 
and heavy carriers and a contains all other contribu- 
tions to the electronic polarizability, e.g., from interband 
transitions; a L is the lattice contribution. Electronic po- 
larizabilities are calculated in the free electron approxi- 
mation: 
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where and are the free-electron energy and Fermi 
function for the i — l,h particles, respectively. 

Neglecting a L for the moment, an approximate expres- 
sion for the dispersion of the APs can be obtained in the 
long- wavelength limit. Using., the first terms of the ex- 
pansions given by LindhardJ£3 
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Vi is the Fermi velocity of the zth particle, and ui^ = 
47re 2 7ii/m* is the plasmon frequency, with m* the effec- 
tive mass and nj the particle density; the imaginary parts 
are neglected. Substituting these in Eq. (uh and defining 



£0 gives for the acoustical plasmon frequency, 
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where v p — Vh y/Nh/3Ni is the acoustical plasmon sound 
velocity, JVf is the density of states of the ith particle 

at the Fermi level, and A = 



— is the screening 
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length of light particles. The Landau damping of acous- 
tical plasmons due to the continuum of light particles is 
neglected. 

In a more complete calculation, we use Eqs. ([]]) and the 
full Lindhard expressions for a 1 . To include the phonon 
contribution we use a simple, approximate expression for 
the lattice polarizability a L 



^M = -^( g ) 2 ( E ° + a q M)/o; 2 , 
and take in the r to A direction 
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where q/qo is the fractional distance to the BZ boundary 
and Co is adjusted to give the frequency of the LA phonon 
in MgB2 at the zone boundary, as taken from Refs. 
||. This implies that the "bare" phonon frequency uj^ h 
already contains the screening due to e° and the light 
carriers, hence this screening is taken out in Eq. (|^). 

The results of such a calculation are shown in Fig. 
||a. The unperturbed or bare phonon dispersion curve is 
given by the dashed line and the acoustical plasmon by 
the dotted solid curve. The latter is Landau damped as 
it enters the heavy particle continuum at q c ~ 0.25. Note 
again that the AP mode shown here is calculated with- 
out taking into account the phonon contribution. Also 
the bare phonon is before the interaction with heavy 
particles is included. The result of the full calculation 
(e = Eq + a 1 + a h + a L ) are shown by the solid line. Af- 
ter the phonon and heavy particles are allowed to inter- 
act, the phonon frequency rises above the unperturbed 
curve and continues to be affected even after the plas- 
mon is damped, causing the anomalous behavior to de- 
velop. Interestingly, the perturbed phonon remains well- 
defined until it enters the heavy particle continuum where 
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FIG. 2: (a) An example of the calculation of the anomaly 
given by the acoustical plasmon mechanism, as described in 
the text, (b) The anomaly as a function of the mass of the 
heavy particle. 



its width increases rapidly in the region of the anomaly 
due to the continuum of heavy particle contribution to 
the imaginary part of e(q,oj). Neutron scattering ex- 
periments should readily pick up this feature if it exists. 
We note that this is also the region where the electron- 
phonon interaction was found to be strong in Ref. |^. 
Some of these effects will be discussed in more detail in 
a subsequent paper. 

To obtain the curves on Fig. |^a we took the mass ratio 
m h/ rn i = 20, and Fermi momentum ratio fcpy/feir/i = 
2.4. ni = 10 21 cm" 3 , n h = 0.75 x 10 20 cm'" 3 , vi = 
2.4 x 10 7 cm/s, and Vh = 5 x 10 5 cm/s. This corresponds 
to the choice of light mass m,* = 1.5too- Also a value 
of £o = 30 was assigned. As was the case with Nb, it 
is difficult to reconcile these parameter values with the 
published band structure calculations and this remains 
a serious drawback for the application of the acoustical 
plasmon concept to lattice dynamics. 

Additional dispersion curves, obtained by changing the 



effective mass of the heavy particles over a limited range, 
are shown on Fig. |b. By varying other parameters, a 
wide variety of behavior can be obtained, which may be 
of interest for other materials. 

We followed Frohlich and Rothwarf a bit further by cal- 
culating T c from a modified BCS-like expression, kT c = 
huipm exp(— 1/F). uj pm was taken to be the frequency 
at which the plasmon enters the heavy hole continuum 
and F was calculated as described in Refs. [I(| , |l7| . Val- 
ues of T c quite close to the experimental value could be 
found but they were very sensitive to the relevant pa- 
rameters which could be changed over a fairly wide range 
without greatly changing the dispersion curves. Having 
in mind the kn own , restrictions on the acoustical plas- 
mon mechanismllall3, we do not consider these results of 
direct, fundamental significance for the high transition 
temperatures in MgB 2 , that is to say, we do not argue 
that superconductivity is produced by acoustical plas- 
mons. However, the results do suggest that the electronic 
structure of this material leads to some resonance-like in- 
teraction of the electrons or holes with "bare" phonons. 
How this might come about from another viewpoint is 
illustrated by the following calculation. 

In Eq. (|Tj), a h , the polarizability of the heavy carriers, 
which gives rise to the acoustical plasmons, was replaced 
by a q-dependent polarization term, a lh , ascribed here 
to interband transitions that were previously included 
in the constant e°. a lh (q) was chosen to be a "skewed 
Lorentzian" of the form 
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We initially fixed go at 0.7 to locate the maximum in 
the vicinity of the dip in the phonon dispersion curve. 
Modification of the bare phonon dispersion and trial and 
error with the parameters a and b readily provided a sat- 
isfactory representation of the anomalous behavior. The 
results of three such calculations for different values of go 
are shown by the upper curves on Fig. |^. We note that 
the importance of interband transitions has been empha- 
sized in recent electronic structure calculations and par- 
ticularly those of Refs. pc|p4"[ 

An interesting aspect of this calculation is that the 
g-dependence of a lh correlates closely with that of the 
"supecshell frequency", M^ s (q), in Weber's double-shell 
modelHy. Smith et alM calculated u ss (q) for TaC by 
assigning the free electron mass to the supershell. The 
supershell motion was then retained explicitly in the dy- 
namical matrix rather than being transformed out as is 
usually done in the shell model. They found that lv ss (g) 
for the longitudinal mode had a pronounced, asymmetri- 
cal dip centered at the position of the anomaly. Loosely 
speaking, the function of the supershell is to provide a 
"resonant polarization" mechanism to screen the relevant 
part of the dynamical matrix. This implies that the shell 
frequency should vary inversely with the square root of 
this polarization. Calculating the g-dependence of u ss 
in this way from Eq. and comparing it with the u) ss 
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FIG. 3: Results for the resonant polarization mechanism as 
described in the text. The upper curves are for the phonon 
and the lower ones are for the supershell frequencies. 



extracted from Ref . |2q gives the results shown in the 
lower part of Fig. |3j; the agreement is obviously quite 
good. The supershell frequencies have been normalized 
to one another at q = and then shifted into the range 
of phonon frequencies. Of course, the introduction of 
a distinct, well-defined mode for the supershell in this 
region of (q, ui) space is a rather drastic simplification 
whose implications are not clear. However, it is remark- 
able that the calculations of Ref. ^4] also point to a sharp 
resonance in the vicinity of 2.5 eV attributed to inter- 
band transitions. In any case, the simple form in Eq. (Q) 
does seem to capture much the same physics contained 
in the double shell model, and perhaps also that in less 
phenomenological treatments. 

The calculations of the phonon spectra of MgB2 
in Ref. which show the anomalies, used a density 
functional theory (DFT), full-potential linear response, 
LMTO approach. The calculations of Ref. |l|, which show 



no anomaly, used a mixed basis pseudopotential method 
and DFTlH. While the results of the two calculations 
are in general agreement, they difiec. strongly in the dis- 
persion of interest here. Kong et aln also calculated the 
electron-phonon interaction and found it to be particu- 
larly large in the vicinity of the anomaly. Since the cal- 
culations of Ref. U obtained the anomalies, their origin 
must be embedded in and, in principle, can be extracted 
from those calculations. Conversely, since the anomalies 
were not obtained in Ref. [IJ a comparison of the two 
approaches is important from both a computational and 
a physical standpoint. As discussed above, attempts to 
identify a well-defined mechanism producing the anoma- 
lies in other materials in the 1970s were inconclusive and 
MgB2 may now provide an impetus to reconsider their 
origins. It will indeed be interesting to see how well these 
calculations fit the experimental lattice dynamical data 
when crystals large enough for neutron scattering become 
available. Of course, any hint of acoustical plasmon ef- 
fects also would be of great interest. 

To summarize, one recent first-principles calculation of 
the phonon dispersion curves in MgB2 shows anomalies 
in the LA branch similar to those found in other high-T c 
electron-phonon BCS superconductors; however, another 
calculation did not give them. These anomalies can be 
calculated using an acoustical plasmon mechanism but 
the band parameters involved seem unrealistic. A res- 
onant polarization mechanism in which interband tran- 
sitions give a q-dependent contribution to the dielectric 
response function is probably more realistic. A detailed 
analysis of the new MgB2 lattice-dynamical calculations 
and comparison with earlier calculations for transition 
metal materials should be invaluable in giving insight 
into the underlying mechanism. 
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